Tunneling of correlated electrons in ultra high magnetic field 
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Effects of the electron-electron interaction on tunneling into a metal in ultra-high magnetic field 
(ultra-quantum limit) are studied. The range of the interaction is found to have a decisive effect 
both on the nature of the field-induced instability of the ground state and on the properties of the 
system at energies above the corresponding gap. For a short-range repulsive interaction, tunneling 
is dominated by the renormalization of the coupling constant, which leads eventually to the charge- 
density wave instability. For a long-range interaction, there exists an intermediate energy range in 
which the conductance obeys a power-law scaling form, similar to that of a ID Luttinger liquid. 
The exponent is magnetic-field dependent, and more surprisingly, may be positive or negative, i. 
e., interactions may either suppress or enhance the tunneling conductance compared to its non- 
interacting value. At energies near the gap, scaling breaks down and tunneling is again dominated 
by the instability, which in this case is an (anisotropic) Wigner crystal instability. 

PACS numbers: 71.10Pm, 72.15Gd, 72.15 



Low-dimensional systems exhibit "zero-bias anoma- 
lies" in tunneling (non-linearities of the current-voltage 
characteristics at small biases), which reflect the renor- 
malization of the density of states by the electron- 
electron interaction. In particular, tunneling into a one- 
dimensional (ID) metal (Luttinger liquid) is character- 
ized by a power-law suppression of the tunneling conduc- 
tance, gx- A three-dimensional metal placed in a strong 
magnetic field that depopulates all but one Landau lev- 
els (ultra-quantum limit, UQL) provides an example of 
a very special quasi- ID system. It is well-knowrJj u that 
repulsive interactions in the UQL lead to a charge-density 
wave (CDW) or a Wigner crystal instability of the ground 
state. This has been confirmed, for example, by exper- 
iments on graphite in high magnetic fieldsH. The most 
complete analysis of this instability for the case of a short- 
range interaction was performed in Ref. [ ^| ] by solving 
the renormalization-group (RG) equation for the inter- 
action vertex. On the other hand, it has recently been 
shown that for the case of a long-range (Coulomb) inter- 
action Qt exhibits a power-law, Lujttinger-liquid-like be- 
havior at sufficiently high energiesQ. This result was ob- 
tained in two ways - via a perturbative Hartree-Fock pro- 
cedure, which results in the RG equation for the trans- 
mission amplitude, and via bosonization in the coher- 
ent state basis - neither of which took into account the 
renormalization of the interaction vertex. In this paper, 
we combine the approaches of Refs. [ |I]j3| ] and [ ^ ] to 
study the behavior of the tunneling conductance in the 
whole energy interval from the Fermi energy down to the 
energy associated with the instability. This is accom- 
plished by solving a coupled system of the RG equations 
for both the interaction vertex and transmission ampli- 
tude. We find that for a long-range interaction there is a 
parametrically wide energy interval in which the flow of 
the transmission amplitude already results in its power- 



law scaling with energy, whereas the renormalization of 
the vertex is not yet important. In this interval, the sys- 
tem behaves as a Luttinger liquid and the results of Ref. 
[ H ] are applicable. At energies close to the gap, the 
power-law scaling breaks down and tunneling is domi- 
nated by the electron properties in the vicinity of the 
critical point. For a short-range interaction, the critical 
point dominates tunneling at all energies, and there is 
no Luttinger-liquid behavior. We also find that above a 
certain magnetic field, when the lowest Landau level is 
strongly depopulated, the interaction enhances the tun- 
neling conductance of the Luttinger-liquid regime above 
its free-electron value. Such an unusual behavior receives 
a natural explanation in terms of electron scattering from 
the Friedel oscillation near the tunneling barrier. 

First we review the procedure of finding the effect of 
the electron-electron interaction on the tunneling con- 
ductance of a metal in the UQLO. Let the magnetic field 
be perpendicular to the contact plane (z = 0) that sepa- 
rates two metallic sides. We consider both the symmet- 
ric configuration, where both sides are in the UQL, and 
the asymmetric one, where one of the sides is made of a 
high-carrier concentration and/or dirty metal so that, in 
the first approximation, the magnetic field does not affect 
that side. The transmission and reflection amplitudes for 
non-interacting electrons, to and tq, are assumed to be 
known. We choose the Landau basis for the free electron 
wavefunctions (non-interacting electrons in the UQL and 
in the presence of the barrier): 

where 

vb Vz (z) = 9(z)t e ip * z + 8(-z)(e lp * z + r e~^ z ) 
yj_ Pz (z) = 9(z)(e-^ z + r e ip » z ) + 9(-z)t e-^ z . (1) 
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(Here and thereafter, we set U = 1 and also define the 
magnetic length as our unit of length: ^jTic/eB = 1). 

Backscattering of magnetically quantized electrons at 
the barrier gives rise to a Friedel oscillation in the charge 
density, whose amplitude decays away from the barrier 
as z~ x . Correspondingly, the exchange, V ex (r, r'), and 
Hartree, Vfj(r), potentials calculated using exact (in the 
presence of the barrier) but otherwise free wavefunctions, 
also exhibit rapid 2fcj?-oscillations and decay as z _1 away 
from the barrier. The first-order correction to the wave- 
function due to interaction is given by 



J dr> J dv"G(r,r';E) [V ex (r',r")+ 
5(r'-r")V H (r")}^ Px (r"), (2) 



where G(r,r';E) is the Green's function in the presence 
of the magnetic field. When the momentum transfer 
(2p z ) matches the wavevector of the Friedel oscillation 
(2k F ), a slow (z _1 ) decay of the oscillations in V ex and 
Vh leads to a logarithmic singularity in S^f(r), and con- 
sequently, in St. For electron's energy close to the Fermi 
energy, we have 



t = t [l-\r \ 2 cT \n(W/E)] 



(3) 



where c = 1(1/2) for symmetric (asymmetric) configura- 
tion, W ~ Ef is the effective bandwidth, To = T(qj_ = 
0), and r(qj_) is expressed via the Fourier transform of 
the interaction potential Uo(q Zl q±): 



r(q±) = 



(2tt) 2 v f 



2vr 

U (2k F ,<U 



- fe i/ 2 c/ (o,k ± ) 
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The interaction potential, Uq, can be both of long- and 
short-range. The first situation corresponds to a single- 
band metal, in which there are no more charge car- 
riers besides those already in the UQL. In this case, 
Uo(q z ,q.i_) = 47re 2 /(g 2 + n 2 ) (for q <C k F ), where K = 
uip/vF and u p is the plasma frequency at B = 0. The 
perturbation theory works if e 2 /vp — k 2 <C 1 , which 
means that the interaction is long range. As the mag- 
netic field reduces the phase-space available for the mo- 
tion along the field, kp and v F decrease with the field (in 
the UQL, k F ,v F oc B^ 1 ). Therefore, the perturbation 
theory breaks down at sufficiently strong field so that 
k ~ 1. The case of a short-range interaction may corre- 
spond to a situation when a small pocket of the Fermi 
surface is in the UQL, whereas other parts of the Fermi 
surface still contain many Landau levels. Screening of 
the interaction among the magnetically quantized carri- 
ers is then mostly due to the "external" carriers, and the 
interaction may be of short-range. 

The first and second terms in Eq.(Q) are the exchange 
and Hartree contributions, respectively. Due to the Pauli 
principle, they enter Eq.(Q) with opposite signs. Note 



that exchange contribution involves integration over the 
transverse momentum , whereas the Hartree contribu- 
tion enters Eq.([|) at q x = 0. This is due to a simple fact 
that the Friedel oscillation in charge density, and hence 
in the Hartree potential, is translationally-invariant along 
the barrier, whereas the exchange potential involves the 
density matrix and is thus non-local. If the interaction 
potential is peaked strongly at q — (e.g., Coulomb 
potential) the Hartree contribution is enhanced due to 
q± = and may in fact dominate over the exchange one. 
In this case r in Eq.(||) is negative and thus the trans- 
misson amplitude is enhanced by the interaction. Such an 
unusual-behavior should be contrasted, e.g., to a strictly 
ID caseQ, in which To = C/o(0) — Uo(2k F ) and is thus pos- 
itive for any "realistic" interaction potential, i.e., such 
that Uo{q) is a monotonically decreasing function of q. 
In our case, due to the 3D nature of the problem, the 
exhange and Hartree contributions are not simply given 
by the corresponding forward and backscattering ampli- 
tudes but also involve averaging over the transverse di- 
rection. As a result, Tq can take either sign. 

At the next order, V ex and Vh are re-calculated using 
the corrected wavefunctions and are substituted back into 
Eq.(^). The higher order corrections to 5^(r) are higher 
powers of logs. Similarly to tunneling, through a barrier 
of weakly interacting electrons in 1DJ2I summation of the 
most divergent corrections to t in all orders of the per- 
turbation theory can be performed via an RG equation 



dt 



-cT t(l-\t\ 2 ), 



(5) 



where £ = In (W/E) and i(£ = 0) = t . The solution of 
(B) is 



t (E/W) 



y/\r \ 2 + \t \ 2 (E/W) 



cr„ 



(6) 



For small values of |io| (high barrier), Eq.(g) reduces to 
a power-law form t = to (E /W) cT ° . The tunneling con- 
ductance is related to the transmission amplitude via the 
Landauer formula. 

The power-law solution for t is obtained using the bare 
interaction vertex. However, the vertex is also subject to 
renormalization. The flow of r(qjr,£) is described by 
the integro-differential RG equationO, which has a more 
transparent meaning when written for the Fourier trans- 
form of 7(r_L,f) = 1/(2tt) 2 / d 2 qi _r(q ± ,C)e- lcl± - r± : 



dt 



x 1 1 — e lr±Ar '^- ) , 



(J) 



where aAb = a x b y — a y b x . The initial condition for Eq. 
(|) is given by r(q x ,£ = 0) = r(q ± ), where r(qj_) is 
defined by Eq. (Q). The flow of the vertex affects the RG 
equation for t, in which now the bare vertex has to be 
replaced by the renormalized one: 
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! = - c r( q± = o,£H (i H*l 2 ) 



(8) 



Eqs. (Q), (|7|) and (g) provide a full description of tunnel- 
ing into a three-dimensional metal in the UQL. 

Few words about the general features of Eq.(Q) are 
now in place. Physically, 7(1" j_, £) plays the role of an in- 
teraction potential between two electrons whose guiding 
centers are separated by a "distance" |r j_ | across the mag- 
netic field. (Note however that 7(1"^, 0) is related to the 
Fourier transform of the bare interaction, cf. Eq. (|]).) 
The first and second terms in the brackets in Eq. (m) cor- 
respond to the Peierls and Cooper interaction channels, 
respectively. A strictly ID case is recovered by putting all 
guiding centers onto the same line: r± = r' ± . In this case, 
r± A r' ± = and the Peierls and Cooper channels cancel 
each other exactly. This cancellation reflects the well- 
known fact that there are no instabilties for ID (spinless) 
electrons and the system remains in the massless phase 
(Luttinger liquid) down to the lowest energies. In the 
UQL case rj_ ^ r' ± , Peierls and Cooper channels do not 
cancel each other, and an instability becomes possible. 

For a short-range interaction, the (^-dependence of Uq 
can be neglected: £/o(0,qj_) — » uq, Uo(2kp,qj_) — > u 2 k F 
andT(q ± ) = (u - u 2 k F ) /(2tt) 2 v f e~«i /2 . For repulsk 
the solution of Eq. (Q) with the Peierls term alone isElt 



r( q± ,0 = (r- 1 (q i )-e)" 



(9) 



The vertex diverges most rapidly at qj_ = and the 
equation for the pole £ = r _1 (0) yields the gap A = 
Wexp(— 1/To). As the electron density remains homo- 
geneous across the field, the instability correspond to 
a CDW along the field. Inclusion of the Cooper term 
merely changes the magnitude of the gaptl. The RG equa- 
tion for [Eq.(§)] with r(q X) £) from Eq.(g) allows for 
an analytic solution: 



t = 



taf(E) 



V\ro\ 2 + \to\ 2 P(E) 



(10) 



where f(E) = [T Q 6{E - A)ln£/A] c . Near the gap 
(E w A), the tunneling conductance gp oc \t\ 2 assumes a 
universal form 



g T oc(E- A) 2c 6(E - A). 



(11) 



Comparing t from Eq.(|Io|) with Eq. (||), obtained with- 
out taking into account the flow of the interaction vertex 
r, we see that the two coincide only at the lowest order 
in the interaction (first log-correction to to). In contrast 
to Eq.(||), the full solution shows that the higher order 
logarithmic terms do not sum up into a power-law. This 
is already seen from the fact that the first log-correction 
to t [Eq.@] becomes of order to at E ~ A, i.e., when 
the renormalization of T cannot be neglected. In other 
words, the CDW critical point dominates the physical 
properties in the entire energy interval from W to A. 



We now turn to the case of a long-range interaction. It 
is instructive to consider first an example of an ultra-long 
interaction: Uq(t) = const — > 7(rj_,0) oc S(r±). To ana- 
lyze the flow of 7(rj_, £) for small £ it suffices to substitute 
7(rj_,0) into the right-hand side of Eq. (]?]), upon which 
the latter vanishes. It means that the ultra-long range in- 
teraction is equivalent to a strictly ID case, when j(r±, £) 
is invariant under renormalization. The example just 
given is certainly not realistic: even for the Coulomb in- 
teraction the ground state is unstable with regard to the 
CDW formation!! However, remnants of this behavior 
survive for the Coulomb potential in a sense that 7(rx, 
flows slower than f (£). As a result, there exists an inter- 
mediate energy scale, Epl, (A <C Epp -C W) such that 
at E ~ Epl the higher order logarithmic terms gener- 
ated by Eq.(|^) have already summed up into a power-law 
form but the renormalization of Y is not significant yet. 
Epl may be defined by the condition that the lowest- 
order correction to to becomes equal to to- 



PL 



hx{W/E PL ) = 1/cIV 



For a screened Coulomb potential, Eq.l 
the condition k <C kp < 1 yields 

r = (e 2 /TTv F ) \\nK\-(4k F 2 ) 



together with 



(12) 



A simple estimate for Epl can be obtained if 
the first (exchange) term in Eq.([lJ) is larger than 
the second (Hartree) one, in which case Epl ~ 
W exp(—Trvp /ce 2 \\n k\). At the same time, for kp ~ 1, 
the gap is estimated as A = W exp(— avp/e 2 ), where 
a ~ 1. We see that £pi > A due to the presence of a 
large logarithm under the exponential. 

As it has already been pointed out, To is not necessar- 
ily positive: if kp is sufficiently small (the lowest Landau 
level is strongly depleted), the Hartree (backscattering) 
term dominates and To < 0. The dependence of To on kp 
is shown in the inset of Fig. [I] for a range of k. Recalling 
that both k F and k are magnetic-field dependent, one can 
express the field Bo at which To changes sign in terms of 
the nominal field Bq , at which the UQL is achieved, and 
the familiar gas parameter, r s , as B ~ 2Bq/ In 1 / 3 CY^ 1 , 
where C ~ 1. For r s ~ 1 Bo ~ 2Bq. The gap as a func- 
tion of kp and n was obtained by solving the RG equa- 
tion (0) numerically with Eq.(|J) as an initial condition. 
The ratio £c/£pl = (\nW/A)/\n(W/E PL ) is plotted in 
Fig. [l[ All numerical results presented here are for a sym- 
metric contact, i. e., c = 1. For a sufficiently long-range 
interaction (cf. filled circles for k = 0.1 in Fig. [l]) there 
is a wide range of kp (and thus of the magnetic field) in 
which £, c /£,pL > 1. Once this condition is satisfied, there 
exists an energy interval A -c E % Epl (Cpl ^ £ ^ £c) 
in which the system effectively behaves as a ID Luttinger 
liquid, exhibiting a characteristic power-law scaling of 
the tunneling conductance. (For Epl <C E <C W the 
power-law reduces to the first log-correction.) Near the 
gap (E ~ A), the power-law behavior crosses over to 
the threshold behavior [Eq. If the interaction is 
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not sufficiently long-ranged, there is no energy interval 
in which £ c /£pl > 1 is satisfied. This is illustrated by 
the open circles in Fig. |l| (« = 0.6). (We did not plot the 
points corresponding to kp < K, as it would have been 
incosistent with our weak-coupling approximation.) 




FIG. 1. Ratio £ c /£pl for k = 0.1 (filled circles) and k = 0.6 
(open circles). The dashed line corresponds to the value 1. 
Inset: To/(e 2 /tyvf) as function of kp for various values of k 
(from top to bottom: k = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6). 
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FIG. 2. Solid lines: solution of the full system of RG 
equations Eqs. (8, 7, 4) for the screened Coulomb potential. 
Dashed lines: solutions of Eq. (8) with To given by its bare 
value. In the main plot to = 0.1, fop = 0.8 and the values of 
k are indicated in the figure. Inset: kp — 0.4. The curves, 
from bottom to top, correspond to k = 0.1, 0.2, 0.3 and 0.4. 



A full solution of Eqs. (Q,0,||) illustrating various 
regimes of the i(i?)-dependence for a screened Coulomb 
potential is represented by solid lines in Fig. || for t = 
0.1. Dashed lines are solutions of Eq. (||) only with Tq 
given by its bare value [Eq. 0], i. e., without taking into 
account the renormalization of the coupling constant. In 
the main panel, kp = 0.8. The region of small £ shows 



a good agreement between the full solution (solid lines) 
and the Luttinger-liquid solution (dashed lines). As E 
approaches A (£ — > £ c ), scaling breaks down and t(E) 
vanishes at E = A. In the inset of Fig. 0, kp = 0.4. The 
rise of t with £ corresponds to the negative value of the 
bare vertex, Tq. This corresponds to an enhancement of 
gr above its non-interacting value. Had the vertex renor- 
malization been absent, t would have reached unity in the 
limit E — > (£ — > oc). However, the vertex renormal- 
ization curbs the rise in t and eventually brings it down 
to zero at E = A. At the instability, the system be- 
comes an insulator. For long-range interaction, we find 
that the vertex diverges most rapidly at finite q±_ (for 
n ^j— jO-j, q± — »w 0.7). This is consistent with previous 
workau and suggests that in this case the instability is 
of the Wigner-crystal type, with charge modulation both 
along the field and in the transverse directions. 

In conclusion, we found that the range of the electron- 
electron interaction plays a crucial role both in the nature 
of a field-induced instability and in tunneling at energies 
above the corresponding gap. For a short-range interac- 
tion, there is a CDW instability which dominates tunnel- 
ing in the whole energy interval. For a long-range inter- 
action, depending on the depletion of the lowest Landau 
level, tunneling may exhibit two types of characteristic 
behavior. For not too strong depletion, the tunneling 
conductance decreases as a power-law of the relevant en- 
ergy scale at higher energies and vanishes in a threshold- 
like manner near the gap (which for the long-range case 
is the Wigner-crystal gap) . For a stronger depletion, the 
conductance first increases as a power of energy, then 
goes through a maximum and finally vanishes at the gap. 
These predictions are amenable to a direct experimental 
verification by tunneling into low-carrier-density materi- 
als, e.g., graphite. 
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